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Abstract 

First order phase transitions in general proceed via nucleation of bubbles. A 
theoretical basis for the calculation of the nucleation rate is given by the homo- 
geneous nucleation theory of Langer and its field theoretical version of Callan and 
Coleman. We have calculated the nucleation rate beyond the thin wall approxima- 
tion by expanding the bubble solution and the fluctuation determinant in powers 
of the asymmetry parameter. The result is expressed in terms of physical model 
parameters. 
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1 Introduction 

First order phase transitions are a common phenomenon in statistical mechanics and in 
field theory [0] . They are characterized by the discontinuous change of an order parameter 
or other physical quantities as some driving parameter, e.g. temperature, is varied. In 
general they are associated with a latent heat. In the theory of elementary particles 
different phase transitions, which play a role in the evolution of the early universe, are 
predicted to be of first order. Among them is the electroweak phase transition, which has 
*present address: SAP AG, D-69190 Walldorf 



been investigated intensively in recent years, and the grand unification phase transition, 
which might be related to the inflationary epoch of the universe [Q. 

There are different mechanisms by which a flrst order phase transition can take place. 
In many cases it proceeds via nucleation of bubbles. Consider evaporation for example. If 
the temperature crosses the transition point the system enters a metastable state. In this 
state bubbles of the new, stable phase form spontaneously, which may then expand and 
lead to the completion of the transition to the new phase. Such metastable states have 
flrst been mentioned by Fahrenheit A theory of the formation of bubbles in liquid 
systems has been developed by Becker and Doring 0] . In the framework of the Ginzburg- 
Landau theory of phase transitions a phenomenological treatment was given by Cahn and 
Hilliard . The theory of bubble nucleation was put on a profound theoretical basis by 
Langer |||, |^, . His approach allows for a systematical treatment of the nucleation rate. 
A review is given in . In the context of quantum fleld theory the nucleation theory was 
developed by Voloshin et al. |]ID[, Callan and Coleman [TT|, and Affleck |]T3[. Callan 
and Coleman presented an approach to the decay of an unstable vacuum in the framework 
of Euclidean quantum fleld theory. Although developed independently, it is very similar 
to hanger's formulation in terms of functional integrals. A nice exposition is given in 



Coleman's book [14 



The starting point of nucleation theory is the classical Ginzburg-Landau potential for 
the order parameter. It has an absolute minimum, corresponding to the stable phase, 
and one (or more) other minima. When the flrst order phase transition is approached, 
a previously higher minimum gets lower and lower, and when the transition point is 
crossed, it becomes the new absolute minimum. There is a barrier between the minima 
such that the system does not immediately go over into the new minimum but remains 
in a metastable state. This state is stable against small fluctuations. Due to fluctuations 
small regions (bubbles) of the stable phase may form spontaneously. Their creation leads 
to a gain in energy proportional to the volume, 

Att 

-nv = -R'v, (1) 

where R is the radius of the bubble and t] is the difference of the potentials between the 
two minima. On the other hand a surface energy 

Hs = ^irR^a (2) 

has to be supplied, and the total energy associated with the bubble is approximately given 
by 

HkiR) =ns + nv = ^7TR\a - ^r]) . (3) 

For small R this function increases with R so that small bubbles tend to shrink back to 
zero. Only if the radius exceeds the critical size of 

Rc = —. (4) 
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where 

<i?M^£)=0, (5) 

the bubble will expand and lead to the transformation of the metastable phase into the 
stable one in the whole volume. 

The process described above is called homogeneous nucleation in contrast to heteroge- 
neous nucleation, where impurities, like dust or ice crystals, trigger the phase transition. 
The cosmological phase transitions mentioned earlier are homogeneous. The formation of 
bubbles in homogeneous nucleation theory is analogous to quantum mechanical tunneling 
through a potential barrier. In fact, the description of tunneling by means of Euclidean 
functional integrals leads to a completely equivalent formalism. This fact is the basis of 
the relation between hanger's work and that of Callan and Coleman. 

The rate in which the phase transition proceeds is essentially determined by the 
average time until a critical bubble forms spontaneously. A critical bubble of radius Rc is 
a solution of the field equations coming from the Ginzburg-Landau Hamiltonian. It has 
an energy Tic = TibiRc)- The nucleation rate F per time and volume is proportional to 
the Boltzmann factor of a critical bubble and can be written as 

r = Ae-^^ , (6) 



which has already been found by Arrhenius [15 



For practical applications it is important to know the prefactor A. hanger's theory 
gives an expression for A in terms of the determinant of the operator of fluctuations 
around the critical bubble. It is the main object of this article to calculate the nucleation 
rate including the prefactor in the framework of scalar field theory, i.e., Ginzburg-Landau 
theory with fluctuations. Some elements of the calculation have been supplied by Langer 
0, but a complete analytical calculation has been missing in the literature so far. A 
numerical method for the evaluation of the nucleation rate has been presented by Baacke 
and Kiselev 

In general it is not possible to find an analytical solution of the field equations for 
finite potential differences rj. An approximation, where the field equations can be solved 
exactly, is the "thin wall approximation" [|1^]. This is the limiting case where rj is much 



smaller than the height of the potential barrier. In this case the thickness of the wall of a 
critical bubble is much smaller than its radius, and its density profile can be approximated 
by a step function. 

In this article we calculate the nucleation rate analytically beyond the thin wall ap- 
proximation. We do this by expanding all quantities in powers of rj and calculating the 
logarithm of the functional determinant of the fluctuation operator in terms of powers 
of rj. The leading term corresponds to the thin wall approximation. For the calculation 
of the determinant we employ the Seeley expansion of the associated heat kernel on the 
one hand, and the spectrum of the fluctuation operator on the other hand. Ultraviolet 



divergencies require renormalization as usual. The resulting expression for the nucleation 
rate F is expressed in terms of renormalized parameters of the effective potential. 



2 Nucleation theory 

As shown by Langer the nucleation rate F, that is the decay probability per time 

and per volume of a metastable state represented by a local minimum of a potential, is 
proportional to the imaginary part of a certain energy: 

F = -2 ImE . (7) 

The energy E is given by the logarithm of a functional integral 

A/" / [d<P] e~^^*^ , (8) 



with appropriate boundary conditions, where (f){x) is the local order parameter and 7i 
is the Ginzburg-Landau Hamiltonian (a factor /c^T has been absorbed into Ti). The 
Hamiltonian is given by 



n{(f)) = jSx 



(9) 



with an asymmetric potential f/. We consider a potential of the type depicted in Fig. 
ID, with a metastable minimum at = 0+ and a stable minimum at = 0_. Following 
Coleman we call the phase corresponding to the minimum at 0+ the false vacuum and 
the one corresponding to the minimum at 0_ the true vacuum. 




Figure 1: The potential U with false 



and true vacuum 



In a semiclassical approach the desired imaginary part of the functional integral can 
be obtained by means of a saddle point expansion around the classical solution which 
corresponds to the transition from the false to the true vacuum 



11, 12 . This solution 



describes a critical bubble. A bubble centered at the origin is represented by a radial 
symmetric function 0c('"), which depends on r = ^x^x^ only. The boundary condition 
at infinity 



lim (brir) 



(10) 



reflects that there is false vacuum outside the bubble. The presence of the true vacuum 
inside the bubble means that the value of the field at the center is near (h_: 



0,(0) ^ 0_ 

Due to differentiability of (j){x) we must have 



dr 



r=0 



The field equation for the bubble solution is 



'Ill 



(12) 



(13) 



dr'^ ' r dr 

If we interpret r as time and as the coordinate of a particle, then this equation equals 
the equation of motion for a point particle in the reversed potential —U with a time- 
dependent friction term. From the form of the potential it is intuitively clear that there 
is a unique value of 0c(O) near 0_, where the particle starts with zero velocity, then rolls 
down the slope and climbs up the other hill to approach its top 0+ asymptotically as 
time goes to infinity. In fact, Coleman, Glaser and Martin ||18| have proved that such a 
radial-symmetric non-trivial solution exists and that it is the one with smallest energy 
apart from the trivial solution = 0+. The qualitative form of the solution is as in Fig. 
0, which shows a cross-section through the critical bubble. 



Figure 2: Profile of the critical bubble. 

For fields 0(x) near the classical solution 0, the Hamiltonian can be expanded up to 
quadratic terms as 



^[0] = + ^ ld'x{<l){x) - Ux))M{<j){x) - 0e(a;)) + . . . 



(14) 



with the energy of the critical bubble, 



K = ^[0c], (15) 

and the fluctuation operator 

M = -d^ + U"{(j),{x)). (16) 

Due to translation invariance the operator Ai has three zero modes proportional to 
d^j^(j)c{x). Furthermore, there is one negative mode, which is related to the metastabil- 
ity of the false vacuum. Namely, as has been indicated in the introduction, expansion or 
contraction of the critical bubble lowers its energy, which means that the corresponding 
mode belongs to a negative eigenvalue of A^. A proof of the existence of a single negative 



mode under rather general assumptions has been given by Coleman |T9 |. 

From the work of Langer and of Callan and Coleman it follows that the functional 
integral under consideration acquires an imaginary part, which is proportional to 

|A_det'Xr'/'e-^^ (17) 

in the Gaussian approximation, where det' is the determinant without negative and zero 
modes, and A_ is the negative eigenvalue oi M.. If multibubble solutions are taken into 
account in a dilute gas approximation, the final result for the nucleation rate is obtained 
as 



nr\"' 1 det'A^ -^'^ 



ncY 

27r/ 



detA^W 

Here the operator Ai^^^ is the Helmholtz operator defined by 

= -d^ + U"{(j)+) . (19) 

The expression (p!8|) is of the form announced in the introduction. In this work the above 
expression, in particular the functional determinant, will be evaluated by field theoretic 
methods. 



3 The bubble solution 

We consider the standard Ginzburg-Landau potential consisting of a symmetric double- 
well term, 

Us = ^{<P'-v')\ (20) 
and an additional asymmetric term: 

U = Us + ^{4>-v) + Uo. (21) 
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The constant 



+ 



(22) 



8v^g IGv^g'^ 

is chosen such that U{(j)-^) = (cf. Fig. [l|). The parameter rj fixes the asymmetry of the 
potential In particular, the difference between the values of the potential at its minima 
0± is 

t/(0+)- W-) = 2r/ + 0(r/3). 
We define a bare mass parameter m in terms of the symmetric part Us 



(23) 



92 



m 



<f)=V 
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The field equation for radially symmetric fields is 

(fcj) 2 



dr'^ r dr Q ^ ' 2v 

Before we turn to a systematic approach to solving this equation let us consider the 
thin wall approximation. Inspection of the field equation in the light of the mechanical 
analogue mentioned in the introduction shows that for small rj the solution is nearly 
equal to — f inside a sphere of radius R and nearly equal to +v outside. The region 
where (p differs significantly from these values is a thin shell of thickness 6. The thin wall 
approximation amounts to 



V 



0, 



(24) 



(25) 



— V 

+v 



for r <R- 9/2 

for R- e/2 <r < R + e/2 

for r> R- 9/2 , 



where the "kink" 



is a solution of 



4>k{f) = f tanh { —{r — R) 



m 



(26) 



(27) 



dr"^ 6 



and 6 = A/m. 

The energy of a bubble solution can be written as 



4-71 drr 



2711] 



drr'^{(j) — v) 



(29) 



The second term is substantially different from zero only inside the bubble. It yields the 
volume contribution 

nv = -—RS- (30) 

The first term gets a substantial contribution only inside the wall, 



An 



R+e/2 
R-9/2 



dr r 



dr 



+ UM + U0 
7 



AttR^ dr 



dr 



AnR^a 



(31) 



with 



9 

m 



a = 2 — . (32) 
9 

For small asymmetries t] the critical radius Rc = 2a /t], for which the energy is stationary, 
gets large. Therefore the wall is indeed thin compared to the size of the bubble. The total 
energy in this approximation is 



— ■ (33) 



For finite 1] the solution of Eq. (|25|) cannot be written in closed form. In our approach 
the solution is constructed by means of an expansion in powers of rj. It is convenient to 
introduce dimensionless variables 

r = ^r, R= —R, ^ = r - R, fj = ^77, (34) 

(^(0 = U{r) , (35) 
where the value of R will be fixed later. The field equation in these variables is 

Based on the thin wall approximation we write a Laurent series as an ansatz for the 
critical radius, 

R = + ao + aifj + a2fj^ + . . . , (37) 
and expand the field equation into powers of fj: 



(Tip 2 ^d(p ^2 2(^ + oo) rfv' 0^2 4 



,,2 + r ' + 2¥^(¥^^ - 1) + T^r^ + 0{n') = . (38) 

at, a_i a4 ci-i 

Its solution is obtained perturbatively up to second order by means of the expansion 

(^ = (^o + ??<^i + r/V2 + 0(f). (39) 
To zeroeth order we get the well-known kink, 

V9o(0 = tanh(0 . (40) 



The field equation to first order fixes the leading coefficient in R, Eq. (|37|) , as 



a_i = l. (41) 
The first order solution obeying the correct boundary condition at ^ — »• cxo reads 

991 = -i - csech^C (42) 
3 
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with a free parameter c. The constant term reflects the shift of the minimum 

^± = ±1-1^1-^ + 0(7?^)- (43) 

The term proportional to c can be traded against a shift in the critical radius R in the 
lowest order solution according to 

tanh(^ - cfj) = tanh ^ - cfj sech^^ - 77^ sech^^ sinh C + O (77^) . (44) 

We can therefore set c = and remain with 

= (45) 

The equation to second order implies 

ao = (46) 

and has the solution 

MO = -|(tanh^-l) + |(coshe-sinhO'-^e'sech^^-^esech^e 

-ln(l + e"^^) [-Csech^^+ itanh^") - i ln(l + e"^«) sinh^ cosh{ 
__Ltanhe + ^sech^e7^(0, (47) 

where we define 

- f rtanh«'- (48) 
Jo 

Whereas the first order solution only corresponds to shifts of the minimum and of the 
critical radius, the second order solution describes true deformations of the bubble. The 
boundary condition at r = 0, i.e. ^ = —R, is fulfilled order by order in fj. For example, 
the leading order solution yields 

^',{-R)=e-'/^{A + Om, (49) 

which vanishes to all orders in 77. Similar observations hold in higher orders. 

With the expression for (p we can calculate the energy of a bubble, which in dimen- 
sionless quantities is given by 



= ^ fdrr^ {{^'iOr + [i^'iO - I? - ivl - If] + [Iv - ^-.)] ] ■ 

(50) 

The integrands are centered around the critical radius f — R. The integration range in 
^ can be extended to the whole real axis. The error coming from this is proportional to 
factors of the type e-'^°^^^-/v and vanishes to all orders in fj. 



9 



Prom the parity of the functions (fikiO it follows that the expression for the energy is 
of the form 

127rm 



He 



9 



-{OoR' + Po) + v{LiR + viR'') + fi\02R' + P2) 



+ f{LsR + V3R^) + 0{t). 
An expression for the critical radius R is obtained from the condition 



~dR 



0. 



(51) 



(52) 



Exphcit calculation of the coefficients leads to two more terms in the Laurent series for 
R: 

-?7 + ■ 77" + 0{ff) , (53) 



fj 6b 



so that the bubble is now completely determined to second order. For the energy we get 



He. — 



127rm 
9 



l_ 2(4 - Ott^) 
9^ ^ 



(54) 



4 The heat kernel of M 

Our main task is to calculate the determinant ratio 

det^A^ 
detA^(o) ' 

which is part of the prefactor in the nucleation rate. In dimensionless variables the 
corresponding operators are 



(55) 



M 



and 



92 



+ 6(^^(0-2 



+ 6(^2 _ 2 . 



Substituting the bubble solution yields for the potential 

v{i) = 6^2(0 - 2 = VoiO + WiiO + fv^iO + 0{f) , 

with the coefficients 

VoiO = -6sech'e + 4, 

\/i(0 = -4tanhe, 
2 

^2(0 = - +C [4tanh{ + 4sinh{ cosh{- (7 + 61n2)sech^{ sinh{] 

— tanh ^ ln(cosh ^) ^6^ sech^^ + 4 cosh ^ sinh ^ + 6 tanh 
-(l + 61n2)tanh^^-41n2 sinh^^ 
+6sech^^ sinh^ T{^) . 
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(56) 
(57) 

(58) 

(59) 

(60) 



(61) 



For the free operator we find accordingly 



= Q^l-2 = ^0^°) + r/V/°) + r^Vi°) + 0{f) , (62) 

with 

Vi'^ = +4, (63) 

= -4, (64) 

Vr = (65) 

The scahng dimensions of Ai and Ai^'^^ are equal, as we have checked with the help of 
their zeta functions (see below). Therefore in the determinant ratio the four removed 
eigenvalues lead to a supplementary factor: 

det^-^ - M V ^"^^^ (66) 
detA^W \my detM(o) ' ^ ^ 

To get the same number of eigenvalues in the numerator and denominator, we also remove 
the four lowest eigenvalues ujq^^ = V-^-* of the free operator M^^\ which are equal to each 
other: 

det'M , (0) 4 det'M 



detM(o) -'"^^ det'MW ^^^^ 
The prime indicates the omission of the four lowest eigenvalues. The remainder is written 
in Schwinger's proper time representation: 

det'M dt„ ^jA^ _jA^(o), 



'"dSTTM-y. fT.'(e--e--). (68, 

In the integrand we recognize the heat kernels exp(— tM) of the operators M and M^^\ 
respectively. The distribution of the large eigenvalues determines the behaviour of the 
heat kernel for small t, whereas the lowest eigenvalues determine the large-t behaviour. 
The strategy for the calculation of the integral over t is to divide the integration range 
into a small-t part, where the heat kernel is approximated by an asymptotic expansion, 
and a large-t part, where the low lying spectrum is employed . 



Let us consider the small-t region first. For small and positive t an asymptotic ex- 



pansion for the heat kernels, the socalled Seeley expansion, exists |^0[. For the trace of 
the heat kernels in D dimensions it is of the form 



Tr(e-*^ - e-**^"') = (47rt)-^/2 g ^n^^ ^gg^ 



-tM ^-tM^°'> ^ 

n=l 

There are various methods for the calculation of the coefficients O^. Our calculation is 



based on the insertion of a plane wave basis in the manner of [^, ^ . By means of partial 
integrations we managed to express the coefficients in terms of the potential V and the 
Laplacean d^, which for a radial symmetric potential like ours depends only on f p4 
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Inserting the potential V{^), Eq. (^8]), and substituting the expression for R, Eq. 
we obtained with the help of Mathematica |]23[ the result 



1 



(47rt)3/2 



fj2 



112 TT 



t- 



160 TT 



832 TT 
T5" 



11392 TT 
315 



3328 TT 
315 



49664 TT 
24255 



+ (-39.4801 1 + 372.46 - 541.384 + 658.823 - 913.886 + 1225.92 - 



(70) 



From the leading term t^^/^ it can be seen that the integral over t diverges at small t. It 
is well known that this divergence is another disguise of the usual ultraviolet divergencies 
of quantum field theory. They can be treated by means of dimensional regularization. It 
is a peculiarity of the three-dimensional case that the regularized expression is identical 



to the zeta-function regularized one without any additional finite contribution p5[. The 
zeta function of an operator A is in general defined by 



1 roo 

Ca{z) = TtA~' = dt Tr(e-^*) 

i (z) Jo 



(71) 



for sufficiently large Re z, where the integral converges, and continued analytically to the 
rest of the complex plane. The zeta-function regularized determinant is then given by 

d 



In det A 



dz 



CAiz) 



(72) 



2=0 



As we are interested in the ratio of two determinants, and want to exclude zero and 
negative modes, we define 



C'(^) 



z-1 T^„'/„-tM 



dt e-' Tr'fe 



— e 



for Re2; > 1 



(73) 



V{z) Jo 

It can be continued analytically to 2 = by separating the first term in the Seeley 
expansion of Tr(e-*^ - e"**^""), Eqs. (I6|),(0): 



C'(^) 



r(^) Jo 



dt e-^ 
1 



^Tr'(e-*^^ 



e-tA^<«')_e(l-t) ^ 



(47rt)3/2 



tOi 



(4vr)3/2 Y{z) z-l/2' 
Now the ratio of determinants can be expressed as 



In 



det^M 
det'M(o) 



d 

dz 



C'(^) 



(74) 
(75) 

(76) 



z=0 



At this point let us check the relative dimensions of the two operators. If the dimen- 
sions of M and M^^^ differ by some number u this means that 



det' AM 



A" 



det'M 



det' AM(o) det' 
12 



(77) 



Representing the ratio of determinants by the derivative of the zeta-function one obtains 

u = C'(0) . (78) 



In the expression (^4D for the analytically continued zeta function the integral is conver- 
gent. Since l/r{z) vanishes at 2; = 0, we find 

C'(0) = , (79) 

so that the operators have equal dimensions as announced above. 

5 Calculation of the determinant 

As mentioned above, the strategy is to split the t-integration into a small-t part and a 
large-t part. Therefore we introduce a parameter A separating the two regions, and split 
the zeta function as 

Ciz) = CA^) + CM (80) 

with 

Ciz) = f'^dt Tr'(e-*^ - e-*^'°') (81) 
T[z) Jo 

and C>(^) correspondingly. In the same way the log of the ratio of determinants is split 



5> 

as 



det'M 

7 = In- 
( 

with 



z=0 



dz 



i> = - -C>{z] 



^3) 
M) 



z=0 

Because the t-integration is not singular for t > A, we can express /> directly as 



I^ = - Tr'(e-**^ - e-**^'°') . (85) 

J A t 



5.1 High-frequency part 

The behaviour of the heat kernel at small t is governed by the high frequencies in the 
spectrum of the operators. Therefore the contribution to the determinant from the inte- 
gration over small t is its high-frequency part. The small-t expansion of the heat kernel 
is given in Eq. (ffOl). For the calculation of C<(^) we have to subtract the contributions of 
the negative mode < and the three zero-modes of M, and to add the contribution 
of the four lowest eigenvalues of M^^^: 
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For Ke z > 1 the integral can be performed: 
1 1 °° n 

C (z) = — = V — A^-y2+n 

^ r(^) {47iy/^ ^,z-3/2 + n 



r(^) Jo ' ' r(z+i) ' Y{z) 

where 

7(a,x) = Te-'r-^rft (88) 





is an incomplete gamma function sec. 6.5]. This expression can be continued analyt- 
ically to 2; = 0, and the derivative at this point yields 

J = 1 ^\ On KU^^nX 

(47r)3/2^^Jn-3/2 / 
+Ei(|cu_|A) - ln(|cj_|A) + 37 + 4Ei(4°^A) + 41n(4°^A) , (89) 

where 7 = — r'(l) = 0.57721 ... is Euler's constant and 

/OO 
—dt , a; > (90) 

Ei(x) = / — cit, a; > (91) 

J X t 

are the exponential integrals |2^, sec. 5.1]. Here the explicit coefficients On from Eq. (|70D 
and the expressions for tu. and lu^'' given in the next section are to be inserted. 

5.2 Low-frequency part 

The low eigenvalues of M and M^°) determine the behaviour of the heat kernels for large 
t. We calculate the low-lying spectrum with the help of a perturbative expansion in fj. 
In spherical coordinates the eigenvalue equation 

Mv{f) = ujv{f), (92) 

goes by means of the usual transformation 

v(r) = - ip(r) (93) 
r 

and the shift 

^ = f-R, = Hr) (94) 

over into 

r ^2 1(1 1 



where V{C,) is given in Eq. (pSf ), and / = 0, 1, 2, ... is the angular quantum number. With 
the help of Eq. (E3) the left hand side is expanded in powers of fj. To lowest order one 
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finds the Poschl-Teller potential Vq = —6 sech^,^ + 4, whose eigenvalues are known exactly 



| 27| . There exist two discrete values 







3 



0, V'o°(e) = \ 7«ech'^' (96) 



V^3°(0 = W-sinhesech^e, (97) 



and a continuum 

ujl = k'^ + A, ken, (98) 
with the corresponding eigenfunctions 

^1! (0 ~ e^^^ (3 tanh^ ^ - 1 - A;^ - 3« A; tanh • (99) 

Because of the radial symmetry, the problem is here only defined along the half-axis f > 0. 
The boundary condition at f = 0, i.e. ^ = —R, is obeyed up to terms which vanish to all 
orders in fj. 

To first order in fj the corrections to the eigenvalues vanish. 

For the second order, we use the following trick [|28|. For every eigenfunction if)^ in 



the zeroeth order there is an operator fulfilling the relation 

[l^n,(-52 + K.)]V'n = ^l^n- (100) 

The second correction in the eigenvalue is then given by 

c.2 = (V;0|\/ifi„ + V2|^°) + /(/ + !). (101) 

This way, we find the second order for the discrete eigenvalues: 

ujI = /(/ + !) -2, (102) 
cul = /(/ + l) + 3-7r2, (103) 



usmg 



Qo = C, ^^3 = 2e-^. (104) 

smh4 



So the discrete eigenvalues of M are 



cuoi = f (/(/ + !) -2) + 0(r/4) (105) 
= 3 + ?72(/(/ + i) + 3_7r2)+0(?7^). (106) 



They are {21 + l)-fold degenerate. 

In particular, the negative mode is given by 



uj_ = ujoo = -2f + 0(77^) , (107) 
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and the three zero modes are cuqi- 

The band of eigenvalues uoi near zero gives a contribution to /> which reads 

dt 



-^(2/ + l)e-"o'*. (108) 

^ 1=2 



This expression can be evaluated with the help of various nontrivial relations involving 
special functions. The result is 

/°(A) = -4t - ^ln(f A) + Co + 0(f A) (109) 

with 

Co = ^ - In 54 - - ^7 = 0.21068 .... (110) 

However, the corrections to the eigenvalues from next order perturbation theory would 
produce corrections of order {fj'^A)^, and therefore we write 

/^(A) = -^-l Hv'^) + oaf A)') . (Ill) 



In a similar way the band of eigenvalues u^i near three gives a contribution 

3 
V 



4(A) = --^r(-l,3A) + 0(r^°). (112) 



with the incomplete gamma function |26, sec. 6.5]) 

T{a,x) = / e-'f-^dt. (113) 



The term of order fj^ has been calculated, but is not displayed, because /° already contains 
an uncertainty of this order. 

Now we turn to the remaining spectrum. The continuous eigenvalues of M and all 
eigenvalues of M*^°^ can be written in the form 

uJki = k^ + V^^^ + fl{l + l) + 0{f,^), k>0, (114) 

with 

1/(0) = 4 _ 4^ _ 1^2 ^ 0(77^) . (115) 
3 

In order to calculate their contribution to the heat kernel one needs the difference of 
the spectral densities gi{k) and gf'\k). We have calculated the spectral densities in the 
framework of perturbation theory. This was done by extracting the phase shifts 6i{k) 
from the asymptotic behaviour of the wavefunctions and using the relation 

It turns out that in n-th order of perturbation theory there are terms proportional to 
(r^i?)", which contribute to the lowest order in fj, because i? ~ 1/fj. Summing up all these 
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terms we have been able to obtain the spectral densities to lowest order only. Omitting 
the details, the contribution of the continuous spectra to the trace of the heat kernels 
(leaving out the four lowest eigenvalues of M'^°^) is given by 

£(2/ + 1) (gi{k) - gf\k))e-^'^' + Ae-'^'"' = (117) 

+ 0(r^°), (118) 
where 

$(x) = ^ /V^'rft (119) 
a/tt jo 

is the error integral. From this we get the contribution to /> by integration over t: 
3r(-l,3A)-^r(-l/2,4A) + ^r(l/2,A)-/^ - $(v/t)e- 

+ 0{ff). (120) 

We now have all contributions to the low-frequency part /> = /° + + /i; available. In 
Fig. ^ they are presented as a function of A for fj = 0.1. One can see that for A > 1 the 
contribution J° from the band near zero dominates the sum. 





5.3 Composition of the determinant 

With all pieces at hand the logarithm of the determinant 

can now be composed. The uncertainty in our result for /< is of order fj'^. On the 
other hand, for /> there are already unknown contributions of order fj^. In the numerical 
evaluations we nevertheless include the constant terms from our calculation of and 
Since the result for / is strongly dominated by /<, as will be seen below, the influence of 
the unknown corrections is expected to be numerically small. 

In an exact calculation I would be independent of the artificial cutoff parameter A. 
Using the approximative expressions above, a dependence on A of course shows up, and 
we have to make a good choice. As a first guide we consider the heat kernel calculated 
from the Seeley expansion, Eq. (^), and from the eigenvalues, respectively. In Fig. § the 
kernel 

i^'(t) = Tr'(e-*^-e-*^^°^) (122) 

is shown as function of t for fj = 0.1 from the two approximations. 

The small-t approximation and the large-t approximation are in good agreement for 
t < 1.3. A value of A near 1.2, where the two curves intersect each other, appears to be 
reasonable for this value of fj. 
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Figure 3: The heat kernel K' as a function of t for f/ = 0.1. The 
dotted curve represents the Seeley expansion and the full curve is the 
approximation from summing over the eigenvalues. 



To make things more quantitative, the various contributions to / and the total sum 
are shown as a function of A in Fig. |. A broad plateau, where the A-dependence is rather 
small can clearly be recognized. 



s o o o 



s o o 



a_ o o o 



- 2_ o o o 



A 



Figure 4: The different contributions to / as a function of A for r] = 0.1: 

/< ( ), /° ( ), J3 ( ), 4( ). The full curve is the sum 

of the four contributions. 

An optimal value for A can be determined in the following way. As discussed above, 
the expression for I can be expanded in powers of fj. The coefficients are functions of A. 
Requiring 

order by order in rj leads to 

A = 1.11672 + 17.3175^2 + 0(7]^). (124) 
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With this value one obtains 



10.158 5 

[V) = 

with an error estimate of 



/(r^) = IZlip. _ r In - 13.63 + 0(77°) (125) 
1]^ 3 



0.0527 



AI{fj) = -- 6.37 + 0(77"). (126) 

The systematic uncertainties of order fj^ come from the corresponding unknown terms in 
/>. Their contribution should be numerically small, as can be seen from Fig. ^. The error 
AJ is estimated from the error of the Seeley expansion, which in turn is taken to be given 
by its highest term. 

5.4 Improved calculation of the determinant 

The result for the determinant can be improved in various ways. First of all, one observes 
that the high frequency (small t) part dominates the result. Also, the low frequency part 
contains uncertainties of constant order in f]. Therefore it is desirable to use as much 
information as possible from the high frequency part, i.e., from the Seeley expansion. 

The determinant could be estimated from the Seeley expansion alone by introducing 
a smooth exponential cutoff [^. With a free cutoff parameter /i one writes 

K'{t) = Tr'(e-*^^ - e"**^'"') = e-'^*Tr'(e-**^ - e"**^*"') e^* , (127) 

and expands the modified kernel as 

Tr'(e-*^ - e-*^'°') e'^* = ^^-^^ + E dnW^ , (128) 

where the coefficients /i„ come from the Seeley expansion ([70|) and the coefficients gn 
from the subtracted negative mode and zero modes. Due to the factor exp(— /xt) the 
t-integration can be extended to infinity. In this way one gets an estimate for I from the 
small-t expansion alone. 

We have used this method and obtained results which are in fair agreement with 
the earlier ones. It is, however, possible to introduce a further improvement. From the 
knowledge of the band of low-lying eigenvalues ujqi it is possible to calculate its contribution 
to / completely The outcome is 



J° = -^Inf + ^-ln54-4C;,(-l) + 0(f) (129) 

= --ln?72 + 1.1727005 + 0(772). (130) 
3 

This result can be employed in the calculation of / by separating the contribution of 
this band of eigenvalues. This means that from the Seeley expansion above the small-t 
expansion of 

E(2/ + l)e--«* = — -- + 0(ft) (131) 
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is subtracted and the rest is treated according to the exponential cutoff method. To the 
result the expression for /° is added finally. 

The value of the cutoff-parameter n has been obtained with the same procedure as in 
the case of A by requiring dl/dfi to vanish. For the log of the determinant we get in this 

way 

Ufj) = - - Inf - 7.98 + Oif) . (132) 

Tj'^ 3 

We have estimated the error by means of the Shanks extrapolation and error determination 
3l| and obtained 

A/(r/) = - 0.084 - 0.795 fi + 0(f) . (133) 

A last improvement is based on the knowledge of the exact leading terms in the small- 
fi expansion of I{f]) |^0[. From the results for the discrete and continuous spectra of M^'^'> 



it is possible to derive 

c 5 

I = 

with 



I = --lln f + 0{fj'), (134) 
rj"^ 3 
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c= — + 31n 3 = 9.9625. (135) 
3 

So we write our final result for the determinants as 

= ^ - "'"-S"' = ^ - 5 1" - 13.52 + 0(,f ) . (136) 

6 The nucleation rate 

Having obtained the determinant of the fluctuation operator the nucleation rate can be 
calculated according to Eq. ([18|) , where the energy of the critical bubble, Eq. ([5^), and 



the negative mode, Eq. (|107|) , have to be inserted. In terms of the mass m and the 
dimensionless parameters fj and 

u = - (137) 

m 

the nucleation rate is 

'327r 1 c 



m 



^3/2 fj7/3 



exp 



U 2 J T]^ 



+ ( ^(Ott^ - 4)^ + 6.845 + 0{u) ) + 0{fj^' 



:i38) 



The parameters appearing in the Hamiltonian are, however, not immediately accessible 
in phenomenological applications or in a field theoretical context. More appropriate are 
the renormalised parameters, which are directly related to measurable quantities. Even 
more important in this context is the fact that in the dimensional regularisation scheme 
around d = 3 dimensions the divergencies in the relation between bare and renormalised 
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parameters are not visible in the one-loop approximation. Therefore physical quantities 
should be expressed in terms of renormalised parameters. 

We shall use the renormalised quantities as, e.g., specified in The renormalised 
mass mn and the field renormalisation constant Z are defined in terms of the inverse 
propagator at small momenta: 



(139) 



The renormalised mass is equal to the inverse of the second moment correlation length 



^(2) ^ ^ _ 
rriR ' 



1 



The renormalised field is given by 



and the renormalised field expectation value 

vn = Z~'/\v + m 



(140) 



(141) 



(142) 



correspondingly. The renormalised coupling, defined in terms of the mass and the field 
expectation value by 

(143) 



9r 



3m% 

2 ' 



has dimensions of a mass. Its dimensionless counterpart is 



9r 
rriR 



The renormalised dimensionless asymmetry parameter is 

Qr 



2m% 



(144) 



(145) 



Up to first order, the relations between the bare and the renormalised quantities are given 
by IH: 

3 



m 



u 



rriR < 1 



-Ur + 0{Uj^) 



1287r 
31 

Ur{1 + Y^'^R + Oiul) 



(146) 



Expressed in terms of the renormalised parameters the nucleation rate is 



m 



R 



3/2 -7/3 
Ur Vr 



exp 



/327r 1 



c 113 



12 



+ 0{ur) 



\ 3 Ur 2 

+ (^{9n' - 4)— + 0.758 + 0{ur)) + 0{fj 

\27 Ur J 



Vr 



(147) 
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This formula is an analytical expression for T which for the first time includes a 
complete treatment of quadratic fluctuations. Compared to the thin wall approximation, 



TWA = exp 



327r 1 



3 URfil_ 



:u8) 



the leading term for small asymmetries rj, 

/32n 1 c 113\ 1 



3/2 -7/3 

Ur Vr 



exp 



\ 3 UR 2 12jfil 



(149) 



completes the thin wall approximation by giving the prefactor in addition to the energy 
of the critical bubble. 

The next-to-leading terms go beyond the thin wall approximation. They can be 
employed to obtain an estimate for the region of validity of the small fjR expansion. 
Comparing the terms proportional to u]^^ in the exponential and requiring the correction 
to be smaller than the leading term we get 

fjR < 0.65 . (150) 

The result for F can be used to obtain estimates for the nucleation rate by inserting 
phenomenological or measured values for the physical parameters. Moreover it can be 
employed in the context of field theory for an estimate of the decay of a false vacuum at 
high temperatures. 

Baacke and Kiselev iTBI have calculated the nucleation rate in three-dimensional scalar 



field theory by evaluating the fluctuation determinant numerically. Our results can, how- 
ever, not be compared directly, because the authors of [T^ pick out the terms relevant in 
a certain high temperature limit. Also the renormalization scheme is different from ours. 

Another numerical calculation of the nucleation rate in the framework of renormalisa- 
tion group improved effective average actions has been presented in [Q. In this case, too, 
it is not possible without further information to compare their results with ours, because 
the results of [|33[ are expressed in terms of parameters, whose relation to the renormalised 



parameters used here is not clear. 

References 

[1] H.E. Stanley, Introduction to Phase Transitions and Critical Phenomena, Clarendon 
Press, Oxford, 1971. 

[2] A.H. Guth, The Inflationary Universe, Addison Wesley, Reading, Mass., 1997 
[3] D.B. Fahrenheit, Phil. Trans. Roy. Soc. London 33 (1724) 78 (in latin). 



[4] R. Becker and W. Doring, Ann. Phys. 24 (1935) 719. 

22 



[5] J.W. Cahn and J.E. Hilliard, J. Chem. Phys. 31 (1959) 688. 

[6] J.S. Langer, Ann. Phys. 41 (1967) 108. 

[7] J.S. Langer, Phys. Rev. Letters 21 (1968) 973. 

[8] J.S. Langer, Ann. Phys. 54 (1969) 258. 

[9] J.D. Gunton, M. San Miguel and P.S. Sahni, in Phase Transition and Critical Phe- 
nomena, Vol. 8, eds. C. Domb, J.L. Lebowitz, Academic Press, London, 1983. 

[10] M.B. Voloshin, LB. Kolzarev and L.B. Okun, Sov. J. Nucl. Phys. 20 (1974) 644, Yad. 
Fiz. 20 (1974) 1229. 

[11] S. Coleman, Phys. Rev. D 15 (1977) 2929, Erratum: Phys. Rev. D 16 (1977) 1248. 
[12] C.G. Callan and S. Coleman, Phys. Rev. D 16 (1977) 1762. 
[13] L Affleck, Phys. Rev. Letters 46 (1981) 388. 

[14] S. Coleman, The Use of Instantons , in "Aspects of Symmetry" , Cambridge University 
Press, 1977. 

[15] S. Arrhenius, Z. Phys. Chem. (Leipzig) 4 (1889) 226. 
[16] J. Baacke and V. Kiselev, Phys. Rev. D 48 (1993) 5648. 
[17] A. D. Linde, Nucl. Phys. B 216 (1983) 421. 

[18] S. Coleman, V. Glaser and A. Martin, Commun. Math. Phys. 58 (1978) 211. 
[19] S. Coleman, Nucl. Phys. B 298 (1988) 178. 

[20] S. Minakshisundaram and A. Plcjcl, Can. J. Math. 1 (1949) 320; 
R.T. Seeley, Proc. Symp. Pure Math. 10 (1971) 288. 

[21] D.L D'yakonov, V.Yu. Petrov and A.V. Yung, Phys. Lett. B 130 (1983) 385, Sov. J. 
Nucl. Phys. 39 (1984) 150. 

[22] R.I. Nepomechie, Phys. Rev. D 31 (1985) 3291. 

[23] Mathematica: A Practical Approach, N. Blachman, Prentice Hall, Englewood Cliffs, 
New Jersey. 

[24] G. Miinster and S. Rotsch, in preparation. 
[25] G. Miinster, Nucl. Phys. B340 (1990) 559. 



23 



[26] M. Abramowitz and LA. Stegun, Handbook of Mathematical Functions, Dover Pub- 
lications, New York, 1972. 

[27] R. Rajaraman, Phys. Rep. 21 C (1975) 227. 

[28] R. Shankar, Principles of Quantum Mechanics, Plenum Press, New York, 1980. 

[29] L. Carson, Phys. Rev. D42 (1990) 2853. 

[30] G. Miinster, in preparation. 

[31] D. Shanks, J. Math, and Physics 34 (1955) 1. 

[32] G. Miinster and J. Heitger, Nucl. Phys. B424 (1994) 582. 

[33] A. Strumia and N. Tetradis, Nucl. Phys. B 542 (1999) 719. 



24 



